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ABSTRACT 

The Very Long Baseline Array has been used at 15 GHz to image the milliarcsecond 
structure of the Seyfert galaxies Mrk 231 and Mrk 348 at two epochs separated by about 
1.7 yr. Both galaxies contain parsec-scale double radio sources whose components 
have brightness temperatures of 10 9 -10 n K, implying that they are generated by 
synchrotron emission. The nuclear components are identified by their strong variability 
between epochs, indicating that the double sources represent apparently one-sided jets. 
Relative component speeds are measured to be ~ 0.1c at separations of 1.1 pc or less 
(for Hq = 65 km s _1 Mpc -1 ), implying that parsec-scale Seyfert jets are intrinsically 
different from those in most powerful radio galaxies and quasars. The lack of observed 
counterjets is most likely due to free-free absorption by torus gas, with an ionized 
density n e £ 2 x 10 5 cm" 3 at T « 8000 K, or n e £ 10 7 cm~ 3 atTw 10 6 6 K, in 
the inner parsec of each galaxy. The lower density is consistent with values found 
from X-ray absorption measurements, while the higher temperature and density are 
consistent with direct radio imaging of the NGC 1068 torus by Gallimore et al. 



Subject headings: galaxies: active — galaxies: individual (Mrk 231=UGC 8058, 
Mrk 348=NGC 262=UGC 499) — galaxies: jets — galaxies: nuclei — galaxies: Seyfert 
- radio continuum: galaxies 



1. Introduction 

Seyfert galaxies have weak and small radio sources, with typical sizes < 500 pc and typical 
centimeter- wavelength powers of <; 10 23 W Hz -1 (e.g., Ulvestad fc Wilson 1989). These sources 
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are apparently produced by weak jets whose axes are determined by the galaxies' central obscuring 



disks, and interacting with thermal gas within ionization cones (e.g., Falcke, Wilson, & Simpson 



1998 



and references therein). The presence of gas disks or tori is supported by the obscuration 
of broad lines in a number of Seyfert 2 galaxies ( [Antonucci &; Miller 1985 ; Miller & Goodrich 
1990| [Kan 1995| ), by VLBI imaging of H 2 masers in NGC 4258 QHerrnstein et al. 1997j ), and by 



parsec-scale HI absorption in other active galaxies ( Peck Taylor 1998 ; Taylor et al. 1999| ) . 



It is of considerable interest to measure the speeds of Seyfert radio jets close to their central 
engines, to attempt to differentiate between "intrinsic" and "environmental" effects. Previous 
measurements of jet speeds in Seyfert cores are rare. In NGC 4151, upper limits of 0.14c and 
0.25c have been measured on scales of 7 and 36 pc ( [Ulvestad et al. 1998] ) . Component positions in 
NGC 1068, measured with the Very Large Array in 1983 QUlvestad, Neff, fc Wilson 1987j ) and with 
the Very Long Baseline Array (VLBA) in 1996 (Roy et al. 1998), imply an upper limit of ~ 0.5c 
for components ~ 20 pc apart. Recently, in III Zw 2, Falcke et al. (1999) inferred an apparent 
speed <j 0.2c on a sub-parsec scale, based on a single VLBI image following a flux outburst. 

The galaxies Mrk 231 and Mrk 348 contain two of the strongest radio sources found in 
Seyferts, and are prime candidates for the measurement of motions in their cores. Within 
larger-scale VLBI structures, each galaxy contains a parsec-scale double radio source (Halkides, 



Ulvestad, fc Roy 1997] ; Ulvestad, Wrobel, & Carilli 1999). Those central sources have now been 
imaged at two epochs using the VLBA;^ this Letter reports measurements of the component 
separation speeds on scales of < 1.1 pc. 



2. VLBA Observations & Data Analysis 

The VLBA QNapier et al. 1994 ) was used to observe Mrk 231 and Mrk 348 at two epochs 
separated by ~ 1.7 yr, between 1996 and 1998. For each galaxy, observing frequencies ranged from 
1.4 to 22 GHz; the images at 15.365 GHz, obtained in left-circular polarization, have the best 
combination of resolution and sensitivity for the compact components, and are the subject of this 
Letter. All data were initially calibrated in AIPS ( van Moorsel, Kemball, fc Greisen 1996| ), then 



iteratively imag ed and self-calibrated in DIFMAP ( |Shepherd 19971 ). O n-source integration times 
ranged from 1.2 to 2.6 hr, resulting in r.m.s. noises of 0.4-1.1 mJy beam -1 for the final uniformly 
weighted images. 

Two-component Gaussian models were fitted in both the (u,v) plane and the image plane; we 
use the image-plane fits here. Flux-density errors (all errors are lcr) were derived by combining 
a 5% scale uncertainty in quadrature with the fitting error (including noise). Estimated size 
errors are 20% in each axis, with point sources taken to have upper limits of half the beam size. 



4 The VLBA is part of the National Radio Astronomy Observatory, a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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Errors in relative component positions are consistent with the measured noises, except for the 
first-epoch image of Mrk 231, where the quoted error is the total range found by using different 
fitting procedures. 



3. Two-Epoch Radio Images 
3.1. Mrk 348 



Mrk 348 has a redshift z = 0.014 ( de Vaucouleurs et al. 1991 ), yielding a scale of 0.31 pc mas 1 
for Hq = 65 km s _1 Mpc -1 (used throughout). This galaxy is a type 2 Seyfert with a hidden 



broad-line region ( Miller &: Goodrich 1990] ), implying that its inner disk is nearly edge-on. It also 



contains a 200-mas triple radio source flNefT fc de Bruyn 1983|) coinciding with optical [O III] 
emission imaged by Capetti et al. (1996). Our VLB A images (Figure ||) show a small-scale 
double source that is aligned with the larger scale radio and optical emission. Phase referencing 
at the second epoch gives a J2000 position for the stronger component of a = 00 ft 48 m 47*1422, 
5 = 31°57 / 25''044, with an error of ~ 12 mas. The relative component separation (see Table [l]) 
increased from 1.46 mas to 1.58 mas in 1.65 yr. Since the stronger component was resolved at the 
second epoch, its centroid could have shifted in absolute position; we estimate a 3a upper limit 
to this shift that is equal to the component size of 0.16 mas. Adding this error in quadrature to 
the nominal position errors, we find a proper motion of 0.073 ± 0.035 mas yr -1 , for an apparent 
speed (in units of c) of /3 app = 0.074 ± 0.035. The implied epoch of zero separation is 19771 20 > so 



the secondary might have had its genesis during a strong flux outburst in early 1982 ( Neff fc de 
Bruyn 1983Q . 



The total flux density of 122 ± 6 mJy at 1997.10 was slightly smaller than the flux density 



of 169 ± 9 mJy measured at 1995.26 ( Barvainis fe Lonsdale 1998j ). Mrk 348 since has undergone 



a major radio flare, with the southern component increasing by a factor of 5.5 between 1997.10 
and 1998.75, strongly suggesting that it is the galaxy nucleus. This is similar to the more extreme 
flares in III Zw 2, which are discussed by Falcke et al. (1999). 



3.2. Mrk 231 



The redshift of Mrk 231 is z = 0.042 flde Vaucouleurs et al. 19"91 ), and the corresponding scale 



is 0.93 pc mas . Mrk 231 is a Seyfert 1/starburst galaxy with a heavily obscured nucleus and 



a total infrared luminosity of ~ 3 x 10 12 Lq ( |Soifer et al. 198*9 ). It contains a 40-pc north-south 
radio source ( Neff fc Ulvcstad 1988 ; Ulvestad et al. 199S| ) embedded within a starburst several 
hundred parsecs in extent ( Bryant &: Scoville 1996 ; Carilli, Wrobel, &: Ulvestad 1998; ). Our VLB A 
images (Figure ^) show a double source with a position angle differing by about 65° from the 
larger scale source. The second-epoch phase-referenced J2000 position (a = 12 fc 56 m 14*2336, 
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5 = 56°52'25''245) is consistent within the 12-mas error with that previously listed by Patnaik et 
al. (1992). The increase in component separation (Table [l]) results in a measured proper motion 
of 0.046 ± 0.017 mas yr -1 , or /3 app = 0.14 ± 0.052, with a zero-separation epoch of 19731 j^. 
Variability by a factor of 2.5 between epochs indicates that the weaker, northeastern component 
is the actual nucleus of the galaxy. 



4. Nature of the One-Sided Sources 



Parsec-scale radio sources have now been imaged in several Seyfert galaxies (e.g., Gallimore 
Baum, & O'Dea 1997; Ulvestad et al. 1998). In NGC 1068, the parsec-scale source is believed 



to represent thermal emission from the accretion torus (Gallimore et al. 1997). However, the 
brightness temperatures in Mrk 348 and Mrk 231 are too high for thermal emission, and instead 
suggest association of the radio components with outflowing jets. Strong variability in one 
component in each galaxy indicates that it is close to the active nucleus, and that the double 
sources represent one-sided jets rather than straddling the galaxy nuclei. 



4.1. Relativistic Boosting? 

One-sided radio structures often are caused by relativistic jets having speed j3c and pointing 



at a small angle 9 with respect to the observer's line of sight (for equations, see Pearson & 



Zensus 1987). We assume an intrinsic spectral index of a = —0.7 (S u cx v +a ) for the off-nuclear 
components. In Mrk 348, therefore, the observed jet /counter jet ratio of R > 17, together with 
Aipp ~ 0.08, can be due to relativistic boosting only if (5 ^ 0.37 and 9 <; 10°. However, the 
inner disk is edge-on ( |Miller Sz Goodrich 199ti| ), and the half-angle of the ionization cone is ~45° 
( Simpson et al. 1996| ); if the radio jet is inside that cone, 9 ^ 45°, inconsistent with motion near 
the line of sight. In Mrk 231, R > 45 and /3 app 0.14, requiring j3 ^ 0.48 and 9 ^ 10° for Doppler 
boosting to account for the one-sided source. If the radio axis is perpendicular to the 100-pc-scale 



disk ( Bryant &: Scoville 1996 ; Carilli et al. 1998), then 9 ~ 45°, also inconsistent with motion at a 



very small viewing angle. 

The above arguments depend on two assumptions: (1) the jets flow along the axes of the inner 
disks, as they apparently do in NGC 1068 (Gallimore et al. 1997) and NGC 4258 ( [Herrnstein et al 



1997| ); and (2) the measured speeds represent the actual jet speeds, rather than quasi-stationary 
structures (such as shocks) through which faster jets flow. The interpretation of fast jets moving 
through slower shocks has been made for Centaurus A, with /3 app ~ 0.1, based on the internal 



evolution of radio components ( [Tingay et al. 1998 ). However, in Mrk 231 and Mrk 348, the 



structures of the off-nuclear components are consistent between the two epochs, so there is no 
similar evidence for higher flow speeds. 
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4.2. Free-Free Absorption 

A straightforward explanation for the one-sided sources is that the "missing" components are 
in the receding jet and are free-free absorbed by ionized gas. The optical depth at frequency v is 

r ff (z/) w 8.235 x 1(T 2 T~ L35 (iz/GHz)- 2 ' 1 E (1) 

( Mczger fc Henderson 196^ ), where T is the temperature in Kelvin and E is the emission measure 



in cm -6 pc. Optical depths of (15 GHz) ^ 4 are needed in order to account for the observed 
jet/counterjet ratios. Assuming a line-of-sight distance of ~ 0.1 pc through the ionized gas, and a 
gas temperature of 8000 K, the average ionized densities would be n e ;> 2 x 10 5 cm -3 at 0.5-1 pc 
from the galaxy nuclei. The resulting column densities of ~ 10 23 cm~ 2 are remarkably consistent 
with the measured X-ray absorption columns in Mrk 348 ( [Smith fc Done 1996 ) and in Mrk 231 



( Nakagawa et al. 1998j ). The density also is consistent with that found for possible absorption in 



a warm, weakly ionized medium in the H20-maser galaxy NGC 2639 ( Wilson et al. 1998j ). On 
the other hand, if the absorption comes from a much hotter gas having T ~ 10 6 ' 6 K, the average 
density required would be n e ^ 10 7 cm -3 . These values are close to those inferred from the image 
of the torus radio emission in NGC 1068 (Gallimore et al. 1997); such disk emission in Mrk 231 
and Mrk 348 is possible, since it would not be detectable by our observations, due to inadequate 
resolution and brightness-temperature sensitivity. 

The free-free-absorption interpretation requires ionized gas densities of 10 5 -10 7 cm~ 3 in the 
inner parsec of Mrk 231 and Mrk 348. These are higher than the values of n e ~ 10 4 cm -3 inferred 
from the absorption ~ 2 pc from the nucleus of 3C 84 ( Levinson, Laor, &: Vcrmculcn 1995| ), and 
n e ~ 10 3 cm -3 inferred for absorption 15-20 pc from the nucleus of Mrk 231 (Ulvestad et al. 1999). 



Thus, our results are consistent with the presence of disks or tori having ionized densities that fall 
gradually from 10 5 -10 7 cm -3 in the inner parsec of the galaxies to ~ 10 3 cm -3 at ~ 20 pc from the 



nuclei. Even though no H2O maser emission has been detected in Mrk 231 or Mrk 348 (Braatz 



Wilson, &: Henkel 1996| ), possibly because of our viewing angle, the one-sided sub-relativistic jets 



are consistent with the presence of megamaser disks. If Tfj(15 GHz) ~ 4, the absorbing gas may 
become optically thin near 30 GHz, and counterjets might be detectable in very sensitive VLBI 
observations at 43 GHz. 



4.3. Comparison to Other Parsec-Scale Sources 



The apparently low jet speeds in Mrk 231 and Mrk 348 are similar to those seen on parsec 
scales in some weak Fanaroff- Riley I ( [Fanaroff and Riley 1974 ) radio galaxies such as Centaurus A 
flTingay et al. 1998j ) and 3C 84 QDhawan, Kellermann, fc Romney 1998 ). However, most such 
objects are relativistic on parsec scales ( piovannini et al. 1998 ), as are most strong radio galaxies 
and quasars ( [Pearson 1996| ). The apparently sub-relativistic Seyfert jet speeds on the same scale 
may be due to interactions with the small-scale gas, or to physics directly related to the energy 
source, such as low black-hole spin rates (e.g., [Rees et al. 1982| ; [Wilson fe Colbert 1995| ). 



- 6 - 



Both Mrk 348 and Mrk 231, as well as NGC 2639 flWilson et al. 1998| ) and III Zw 2 
(Falcke et al. 1999), have radio spectra that peak at 10 GHz or higher. Thus, the Seyfert radio 
sources have general characteristics similar to compact symmetric objects (CSOs), which display 
gigahertz-peaked spectra, symmetric radio sources on scales of ~ 100 pc, and one-sided jets on 
scales of 10 pc ( Taylor, Readhead, &: Pearson 1996 ). Small-scale Seyfert jet speeds also are similar 
to those measured in CSOs at 50-100 pc from their nuclei ( |Owsianik fc Conway 1998| ; |Owsianik 



Conway, &: Polatidis 1998 ). The CSO jets are one-sided at 15 GHz (Taylor et al. 1996), and 



even at 43 GHz (Taylor, private communication), on scales of 10 pc; at least one CSO also shows 
patchy HI absorption on similar scales ( [Peck, Taylor, Sz Conway 1999 ). We speculate that free- free 
absorption with column densities near 10 cm -2 might cause one-sided structures in some CSOs, 
if their disks are larger and denser than in Seyferts. High-energy X-ray spectral studies could be 
used to search for absorption due to this gas. 
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California Institute of Technology, under contract with the National Aeronautics and Space 
Administration. 



- 7- 



REFERENCES 

Antonucci, R. R. J., & Miller, J. S. 1985, ApJ, 297, 621 
Barvainis, R., & Lonsdale, C. 1998, AJ, 115, 885 
Braatz, J. A., Wilson, A. S., & Henkel, C. 1996, ApJS, 106, 51 
Bryant, P. M., & Scoville, N. Z. 1996, ApJ, 457, 678 

Capetti, A., Axon, D. J., Macchetto, F., Sparks, W. B., & Boksenberg, A. 1996, ApJ, 469, 554 
Carilli, C. L., Wrobel, J. M., & Ulvestad, J. S. 1998, AJ, 115, 928 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, R. J., Paturel, G., & Fouque, P. 
1991, Third Reference Catalogue of Bright Galaxies (New York: Springer) 

Dhawan, V., Kellermann, K. L, & Romney, J. D. 1998, ApJ, 498, Llll 

Falcke, H., Wilson, A. S., & Simpson, C. 1998, ApJ, 502, 199 

Falcke, H., et al. 1999, ApJ, 514, L17 

Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P 

Gallimore, J. F., Baum, S. A., & O'Dea, C P. 1997, Nature, 388, 852 

Giovannini, G., Arbizzani, E., Feretti, L., Venturi, T., Cotton, W. D., Lara, L., & Taylor, G. 

B. 1998, in ASP Conf. Ser. 144, IAU Colloquium 164, Radio Emission from Galactic and 
Extragalactic Compact Sources, ed. J. A. Zensus, G. B. Taylor, & J. M. Wrobel (San 
Francisco: ASP), 85 

Halkides, D., Ulvestad, J., & Roy, A. 1997, BAAS, 29, 1375 

Herrnstein, J. R., Moran, J. M., Greenhill, L. J., Diamond, P. J., Miyoshi, M., Nakai, N., & Inoue, 
M. 1997, ApJ, 475, L17 

Levinson, A., Laor, A., & Vermeulen, R. C. 1995, ApJ, 448, 589 

Mezger, P. G., & Henderson, A. P. 1967, ApJ, 147, 471 

Miller, J. S., & Goodrich, R. W. 1990, ApJ, 355, 456 

Nakagawa, T., Kii, T., Fujimoto, R., Miyazaki, T., Inoue, H., Ogasaka, Y., Arnaud, K., & Kawabe, 
R. 1998, in IAU Symposium 186, Galaxy Interactions at Low and High Redshift, ed. D. B. 
Sanders & J. Barnes (Dordrecht: Kluwer), 103 

Napier, P. J., Bagri, D. S., Clark, B. G., Rogers, A. E. E., Romney, J. D., Thompson, A. R., & 
Walker, R. C. 1994, Proc IEEE, 82, 658 



-8- 



Neff, S. G., & de Bruyn, A. G. 1983, A&A, 128, 318 

Neff, S. G., & Ulvestad, J. S. 1988, AJ, 96, 841 

Owsianik, I. & Conway, J. E. 1998, A&A, 337, 690 

Owsianik, L, Conway, J. E., & Polatidis, A. G. 1998, A&A, 336, L370 

Patnaik, A. R., Browne, I. W. A., Wilkinson, P. N., & Wrobel, J. M. 1992, MNRAS, 254, 655 

Pearson, T. J. 1996, in ASP Conf. Ser. 100, Energy Transport in Radio Galaxies and Quasars, ed. 
P. E. Hardee, A. H. Bridle, & J. A. Zensus (San Francisco: ASP), 97 

Pearson, T. J., & Zensus, J. A. 1987, in Superluminal Radio Sources, ed. J. A. Zensus & T. J. 
Pearson (Cambridge: Cambridge University Press), 1 

Peck, A. B., & Taylor, G. B. 1998, ApJ, 502, L23 

Peck, A. B., Taylor, G. B., & Conway, J. E. 1999, ApJ, in press 

Rees, M. J., Phinney, E. S., Begelman, M. C, & Blandford, R. D. 1982, Nature, 295, 17 

Roy, A. L., Colbert, E. J. M., Wilson, A. S., & Ulvestad, J. S. 1998, ApJ, 504, 147 

Shepherd, M. C. 1997, in ASP Conf. Ser. 125, Astronomical Data Analysis Software and Systems 
VI, ed. G. Hunt k H. E. Payne (San Francisco: ASP), 77 

Simpson, C, Mulchaey, J. S., Wilson, A. S., Ward, M. J., & Alonso-Herrero, A. 1996, ApJ, 457, 
L19 

Smith, D. A., & Done, C. 1996, MNRAS, 280, 355 

Soifer, B. T., Boehmer, L., Neugebauer, G., & Sanders, D. B. 1989, AJ, 98, 766 
Taylor, G. B., O'Dea, C. P., Peck, A. B., & Koekemoer, A. M. 1999, ApJ, 512, L27 
Taylor, G. B., Readhead, A. C. S., & Pearson, T. J. 1996, ApJ, 463, 95 
Tingay, S. J., et al. 1998, AJ, 115, 960 
Tran, H. D. 1995, ApJ, 440, 597 

Ulvestad, J. S., Neff, S. C, & Wilson, A. S. 1987, AJ, 93, 22 

Ulvestad, J. S., Roy, A. L., Colbert, E. J. M., & Wilson, A. S. 1998, ApJ, 496, 196 

Ulvestad, J. S., & Wilson, A. S. 1989, ApJ, 343, 659 

Ulvestad, J. S., Wrobel, J. M., & Carilli, C. L. 1999, ApJ, 516, in press 



-9- 



van Moorsel, G., Kemball, A., & Greisen, E. 1996, in ASP Conf Ser. 101, Astronomical Data 

Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes (San Francisco: ASP), 37 

Wilson, A. S., & Colbert, E. J. M. 1995, ApJ, 438, 62 

Wilson, A. S., et al. 1998, ApJ, 505, 587 



This preprint was prepared with the AAS IAT^X macros v4.0. 



-10- 



Table 1. Results of 2-Component Gaussian Fits 



No. 


Flux Density 


Size T 6 a 


Offset 


PA 




(mJy) 


(masxmas, deg) (10 9 K) 


(mas) 


(deg) 






Mrk 348. 1997.10 






1(S) 


96 ±5 


Unresolved > 6.7 ± 1.9 


0.000 




2(N) 


26 ±2 


Unresolved > 1.8 ± 0.5 
Mrk 348, 1998.75 


1.460 ± 0.009 


-16 ±1 


1(S) 


552 ± 28 


0.16 x 0.11, PA -9 ±4 238 ±68 


0.000 




2(N) 


17 ± 1 


Unresolved > 1.8 ±0.5 
Mrk 231, 1996.94 


1.581 ±0.021 


-15 ±2 


1 (NE) 


17 ± 1 


Unresolved > 2.2 ± 0.6 


0.000 




2 (SW) 


51 ±3 


0.32 x 0.18, PA -74 ±4 6.9 ± 2.0 
Mrk 231, 1998.71 


1.081 ±0.030 


-115 ± 1 


1 (NE) 


44 ±3 


Unresolved > 6.0 ±1.7 


0.000 




2 (SW) 


60 ±3 


0.30 x 0.13, PA -88 ±4 12.0 ± 3.4 


1.162 ±0.004 


-112 ± 1 



a All brightness temperatures are expressed in the source rest frames, by multiplying the 
observed brightness temperatures by (1 + z). 
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MRK348 15.365 GHz 




0.5 pc - 



1997.10 1998.75 




2 1 0-1-2 

MAS 

Fig. 1 — VLBA 15-GHz images of Mrk 348 at epochs 1997.10 and 1998.75. The images have been 
rotated by 15° from the cardinal orientation, aligned at the southern radio component, then offset 
from each other horizontally by 2 mas. North and East are labeled, as are the two radio components. 
Logarithmic contours start at 4 mJy beam -1 and increase by factors of 2 to 512 mJy beam -1 . The 
common restoring beam is is 0.80 x 0.45 mas in PA —15°. 
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MRK231 15.365 GHz 
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Fig. 2.— VLBA 15-GHz images of Mrk 231 at epochs 1996.94 and 1998.71. The images have 
been rotated by 25° from the cardinal orientation, aligned at the northeastern radio component, 
then offset from each other vertically by 2 mas. North and East are labeled, as are the two 
radio components. Logarithmic contours start at 2 mJy beam -1 and increase by factors of 2 to 
32 mJy beam -1 . The common restoring beam is 0.55 x 0.42 mas in PA 0°. 



